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Abstract 

Characterization of the American put option price is stih an open issue. From the 
beginning of the nineties there exists a non-closed formula for this price but nontriv- 
ial numerical computations are required to solve it. Strong efforts have been done to 
propose methods more and more computationally efficient but most of them have 
few mathematical ground as to ascertain why these methods work well and how 
important is to consider a good approximation to the boundary or to the smooth 
pasting condition. We perform an extension of the American put price aiming to 
catch weaknesses of the numerical methods given in the literature. 

Key words: Econophysics, American put option, Computational methods, 

Black-Scholes, Option pricing 

PACS: 89.65.Gh, 02.50.Ey, 02.60.-x, 05.10.-a 



In this paper we study approximations for the price of the American put op- 
tion. In recent years there has been an explosion of literature regarding this 
problem. The mathematical literature has focussed on the study of properties 
of the so-called optimal boundaries and the convergence properties of the bi- 
nomial method of pricing (see Kim [1], Jacka [2] and Carr [3]). In the more 
numerical literature, on the other hand, the main objective is to achieve a fast 
calculation of put prices, somewhat disregarding the approximation of the 
optimal boundary and the so-called smooth pasting condition. Huang et al. 
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wrote: "It seems that in order to have a good approximation of the price is not 
necessary to have a good approximation of the optimal boundary" [4]. Subse- 
quent works have improved their methodologies -see for instance Refs. [5,6]-. 
A common feature though among the computational articles is that they never 
investigate why a proposed method work well or not, and how important is to 
accurately approximate the optimal boundary. Aiming to answer these ques- 
tions, our objective is to bridge the gap between these computational investi- 
gations and the more theoretical literature with the hope to bring some new 
insights. We must admit that this is a difficult task and the present study is 
limited to provide few insights. At this step, we solely sketch some results of 
a more extensive work which is still in progress [7]. 

Let us provide the general framework. We consider the classical Black-Scholes 
model consisting of two assets: one risk free bond B whose dynamics is de- 
scribed by 

dB (t) = rB it) dt with B (0) = 1, (1) 

and one stock of price S following the stochastic differential equation ^ 

dS{t) = rS{t)dt + aS{t)dW{t) with S{Q) = Sq, (2) 



and where dW{t) — ^{t)dt is the Wiener process with zero mean and unit 
variance. 

The American put option gives to its owner the right (and not the obligation) 
to sell an asset at a given strike price K at any time up until maturity T. This 
last condition makes the pricing problem very complicated, much more than 
for European options where exercising the option is only possible at maturity 
date T. Under the Black-Scholcs model determined by the pair of Eqs. (1) 
and (2), the American put price at time t with strike price K is defined by 

p^(t, S) ^ sup e-'-(-*) {[K - S(Tr\S(t) = S) , (3) 

where St^T is the set of all stopping times with values in [t,T]. This definition 
involves an expectation over the put payoff, which solely takes the positive 
contribution of the difference between strike K and stock S. We will fix the 
fair price as the highest average of the collection of expectations with stopping 
times r G [t, T]. Hence, the problem reduces to an optimal-stopping-time prob- 
lem and main purpose is to know the best moment to exercise the American 
option. 

^ For sake of simplicity, in this paper we assume no dividend yield although similar 
conclusions can be obtained with non-zero dividend yield [7] . 
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The pricing problem was solved more than a decade ago departing from dis- 
crete time and recursively backward with non arbitrage arguments [1]. The 
so-called binomial pricing method gives a non-closed solution for the Ameri- 
can put option in terms of the optimal exercise boundary c{t) and reads (for 
alternative derivations see Refs. [2,3]) 



T 

p^{t,S)^p^{t,S)+rK J e-'^^''-^^N[-d2{S,c{u),u-t)]du (4) 

t 

where 



S) = Ke-'^'^-'^N[-d2{S, K,T- t)] - SN[-di{S, K,T - t)] (5) 
is the European put option price and 

N{a) e-^'/'dx (6) 

oo 



is the probability integral function whose arguments are 



, ln(x/y) + (r + a'^/2)t , , r 

di{x,y,t) = ^ V d2 = di-aVt. (7) 

ay/t 

The formula (4) says that the price of the American put is the same as the 
European put but with two extra contributions quantifying the premium as- 
sociated to having the right of exercising early the option. In the put case, the 
optimal boundary characterizes the values below which we should exercise the 
American option {i.e., if S{t) < c{t)) and it is solution to the equation: 



T 

K - c{t) = p^[t, c{t)] +rK j e-''("-*)Ar[-d2(c(t), c{u),u - t)]du. (8) 

t 

This integral equation uniquely characterizes the optimal boundary among all 
non-decreasing continuous functions with values in (0, K] [2] and jointly with 
Eq. (4) constitute the non-closed formula for the American put price. Still, 
there is httle hope of finding an analytical and explicit solution to Eq. (8) and 
subsequently to Eq. (4). For practical purposes though it gives us the possibil- 
ity to numerically generate the optimal boundary, from which an approximate 
put price and replicating strategy can be derived. Each method will have their 
own clever tricks in order to make more efficient the computation. 
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Fig. 1. Computation of the optimal exercise boundary in terms of T — t. Eq. (8) 

is taken in its discrete form assuming that the optimal curve follows a stepped 
function and with K = 100. Left hand side plot shows the quick convergence of the 
optimal curve when r = 5% and a = 20%. No big differences exist between the case 
of N = 500 steps and the N = 10, 000 case. Right hand side shows the results for a 
set of r's, and cr's with N = 500 steps. Both graphs take time units in years. 

An alternative method for numerical computation arises from the Black-Scholes 
partial differential equation (PDE) approach, an approach also known as the 
Stefan problem. Taking finite differences, this method asserts that when un- 
derlying is above the optimal exercise boundary [x > c{t)) the PDE reads [8] 

dp^{t,x) a'^x^d^p^{t,x) dp^{t,x) ^ n n (o\ 



while below the boundary {0 < x < c{t)) we immediately exercise the option 
with benefit 

p^{t,x) ^ K -X. (10) 



For continuity and theoretical arguments [8], both price and first order deriva- 
tive over the price (p^ and p^) must be continuous on the boundary curve c{t). 
In particular, last requirement gives the so-called smooth pasting condition. 
That is: 



_ dp^{t,x) 



dx 



-1. 



(11) 



--c{t) 



It is somewhat deceiving that the pricing method mentioned above [1] does not 
clearly show how the smooth pasting condition (11) is incorporated into early 
exercise premium formula (4). This is a serious drawback of the numerical 
methods which are based on premium formula (4). 

Let us study this particular point in more detail. First we insisit that our 
interest is not the efficiency of the computation but the goodness of the price 
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Fig. 2. Computation of the put price derivative on the optimal exercise boundary 
and its resulting error in terms of T — t when K = 100. Left hand side plot shows 
numerical integration (12) with the estimated values of c{t) with a stepped function. 
This estimator of c{t) has non constant derivative and differs from the theoretical 
value -1. Right hand side plot depicts error (15) quantifying the non-fulfillment of 
the smooth pasting condition when S = 95. 

formula (4) based on the smooth pasting condition. We focus on the particular 
case given by Ref. [1] where we use a step function to estimate the theoretical 
optimal exercise boundary c{t). Equation (8) is easily solved by the Newton's 
method due to the monotonicity of the boundary curve. Figure 1 shows the 
rapid convergence of the estimated boundary curve and also depicts the results 
for several values of r and a. Once we get an estimation of the optimal exercise 
boundary, we calculate the American option price given by Eq. (4). 

And, what about the smooth pasting condition? Prom Eq. (4), we can calculate 
the derivative with respect to the stock and evaluated when S > c{t) : 



^ ( ) 

pi{t,S)^p^s{t,S) - ^ I ^Zl- N'[d2{S,d{u),u-t)]du, (12) 

o J a^Ju — t 
t ^ 

where pf(t, 5*) = —N[—di{S, K, T—t)]. This derivative, popularly called delta, 
is a very important quantity since it prescribes the hedging strategy to be fol- 
lowed. In case we deal with the optimal exercise boundary c{t), the smooth 
pasting condition (11) prescribes that this derivative is -1. However, in case 
we take our estimation of c{t) represented by 5" = d{t) and numerically calcu- 
late Eq. (12), we find that our estimation does not obey the smooth pasting 
condition. Figure 2 shows that the numerical values of Eq. (12) are different 
from -1. This questions the validity of the numerical computation and suggest 
a careful investigation (see Ref. [7] ) . 

Having this inconsistency in mind, we provide a more general framework where 
not only optimal curves are allowed. Let us denote by Co,r the set of non- 
decreasing continuous functions d : [0,T] — > {0,K] with d{T) = K. Hence, 
the set Co,T consists of all candidates for the optimal boundary. Consequently, 
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following the spirit of Eq. (3), the price of the American put can be written 



as 



p^{t, S) = sup (e-'^^*''^-'^[K - S{r)]+\S{t) ^ s) = sup p'^{t, S), (13) 



where Tt^d = inf{s > t : S{s) < d{s)} A T. In Ref. [7], we have analyzed the 
time t price p'^{t, S) for a general boundary d e Co,t- The most important 
result found is that determinsitic function p'^ can be expressed as 



p'^it, S) ^p^{t, S)+rK J e-'-("-*)Ar [-d2{S, d{u),u - t)] 



du 



t 



2 J ^^/u — t 



d{u)N' [-d2{S, d{u),u - t)] (l + p^^ {u, d{u)^)) du, 

(14) 



where prime in the third term denotes the derivative dN{x)/dx. In order to 
derive this extension of the early exercising premium formula for an American 
put, wc have required p'^ to be the solution of the PDE (9) with conditions 
Eqs. (10) and (11). Details of this calculation are given in Ref. [7]. The expres- 
sions (9)-(ll) are essential new ingredients and, in some sense, they introduce 
a third term containing the partial derivative over the stock and evaluated on 
the curve d{t). Note that, in case we deal with the optimal exercise bound- 
ary curve c{t), the smooth pasting condition make this term to be zero and 
our extension of the early exercising premium formula would collapse to the 
American option put price (cf. Eq. (4)). Hence, intuitively, the second integral 
in Eq. (14) is able to quantify the error committed in computations aiming to 
get a price for the American put and reads 



Error =~ [ ^ , d{u)N' [-d2{S, d{u),u - t)] (l + pi {u, d{uf)) du. 
^ J \/u — t ^ ' 

(15) 

This quantity is directly related to the smooth pasting condition and lets us 
go beyond in judging the quality of an estimation of the optimal boundary 
based on the number of steps in the time discretization since, as we have seen 
in left hand side plot of Fig. 2, even for a large number of steps the smooth 
pasting condition may not be accomplished. From Eqs. (14)-(15), we thus take 
(Error/p^) x 100 and represent it in the right hand side plot of Fig. 2. Wc have 
done it in terms of the maturity time, for several values of r and a, and when 
S = 95. We observe that error is negative and non-negligible making always 
cheaper the corrected American put. 
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In this paper, we have approached to computational pricing problem using 
tools from mathematical finance. From these investigations, we may assert 
that the important but ignored properties of the American put option distort 
numerical results obtained via the binomial pricing method. We have focussed 
on the simplest numerical technique that takes the optimal exercise boundary 
curve as a step function and showed that the resulting price does not ac- 
complish the smooth pasting condition. We have proposed an extension of the 
early exercise premium formula that captures the effect of the non-fulfilment of 
the smooth pasting condition. We have observed that errors are non-negligible 
and, most importantly, they can be quantified analytically. For this reason, we 
believe that the extended formula (14) can provide fruitful results with a good 
mathematical ground as to ascertain whether a numerical method works well 
and it is thus able to evaluate the quality of existing computational methods 
in the literature. For more results we refer to the investigation in Ref. [7]. 
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